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Abstract 

Using the inversion relation method, we calculate the ground state energy for the lattice 
integrable models, based on baxterization of multicolored generalization of Temperley-Lieb 
algebras. The simplest vertex and IRF models are analyzed and found to have a mass gap. 

1 Introduction 

Looking for new solutions of Yang-Baxter equation is an interesting and important problem in 
exactly solvable statistical models. Some time ago set of new solutions of Yang-Baxter equations 
were obtained [1] by "baxterization" of multicolored generalization of Temperley-Lieb (TL) al- 
gebras - Fuss-Catalan algebras - discovered by Bisch and Jones [2]. These solutions are really 
new and are not equivalent to a fusion of known models. In [1] these solutions were formulated 
in terms of generators of this new algebra, and in principle the Boltzmann weights may be in- 
terpreted both as a vertex model, or as an interaction round the face (IRF) models defined on a 
tensor product of admissibility graphs corresponding to different colors of the TL algebra. Phys- 
ical interpretation of these new integrable models may be different. In addition to construction 
of new integrable vertex and IRF models, an obvious possible realization of these new solutions 
is multicolored generalization of dense loop models [1], polymers and other subjects of exactly 
solvable statistical models (for a review see e.g. [3]). Another new integrable model of the same 
class was obtained as integrable solution for Lorentz lattice gases [4], [5] and dimerised coupled 
spin-1/2 chain [6]. 

In spite of traditional form of the R-matrix, it turns out that the usual algebraic Bethe ansatz 
solution, at least in its naive and straightforward version, reveals some technical problems. One of 
the possible reasons for that might be a non trivial structure of pseudo vacuum, which is hard to 
guess in a way effective for realization of algebraic Bethe ansatz program. The alternative method 
for investigation of physical properties of integrable models, such as ground state energy, its free 
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energy, spectrum and correlation functions is the method of inversion relations. It doesn't require 
explicit form of BA equations or transfer matrix eigenvalues, valid in the thermodynamic limit. 
It was taken up by Baxter (see, e.g. [7], [8]) and further developed by others (see e.g. [9], [10]). 
It exploits functional inversion relations for the transfer matrix in the ground state, and transfer 
matrix analytical properties, in order to calculate ground state energy and low lying excitations 
spectrum [10]. Assumptions on the transfer matrix eigenvalue analytical properties usually should 
be confirmed numerically. The advantage of this method is that it doesn't require knowledge of 
Bethe ansatz equations and eigenvalues of transfer matrix explicitly, but gives a possibility to 
find basic physical quantities of the model quite easily. 

Plan of the paper is the following. We start Sect. 2 with description of the Boltzmann weights 
of the models in terms of multicolored TL generators and their algebraic relations. We analyze 
possible physical regimes, where Boltzmann weights are real and positive, and discuss possible 
vertex and IRF interpretation of the models. Some interesting factorization property of admissi- 
bility graph for IRF models is pointed out. In Sect. 3 we derive functional inversion relations for 
the transfer matrices of two colored vertex model and for the partition function of two colored 
IRF model. We solve them analytically as a product of analytical non zero function, and a func- 
tion with poles or zeros. Excited states are discussed and found to lead to a mass gap. In Sect. 4 
generalization of inversion relations for multicolored models is discussed. We conclude by Sect. 5 
with a discussion of the obtained results. 



2 Colored TL lattice models. 

In this section we recall the results of [1] and demonstrate them for few examples of two, three 
and four colored models. For details of Fuss-Catalan algebras and their representations we refer a 
reader to [2] and [1]. We also discuss vertex and IRF representations of baxterized Fuss-Catalan 
models and find their possible physical regimes. 

2.1 Solutions of Yang-Baxter equation with Fuss-Catalan symmetry. 

New class of TL like algebras found by Bisch and Jones [2] can be considered as a kind of quotient 
of tensor product of few TL algebras. Generators of these algebras may be formulated in terms 
of TL generators defined on a lattice. The relations on generators depend on the parity of the 
site i where the operator acts. The main commutation relations [2] are the following. 

U- '(]{' = Uf'U> = pi(mm(m,pj)U> ,y " (1) 

where 

, n ' otioi2...a i — even 

Pi\P) = \ ■ aj ( 2 ) 

L a k a k -i...a k+ i_ p i - odd 

Here k is the number of colors is, and in the standard strings pictorial interpretation of the 
operators [/( m ) weight a m is attached to each loop of color m. Locality of operators U expressed 
as 

U\ m) uf = Ufu\ m) (3) 

when \i — j\ > 1, or when \i — j\ = 1 and m + p < k. Non commuting operators satisfy cubic 
relations 

TT (m) TT (p) TT (q) n s I U> 'U- ±1 m>q ( 

U i Ui±iUr = Pi(k - p) \ fe _ ra J ?) (4) 
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These are basic relations of the Fuss-Catalan algebra FC k ^ N+ ^(ai, ct2, ■ ■■a k ) which describes 
coils of k(N + 1) strings of k colors. Boundary conditions are fixed by requirement that both 
on the top and on the bottom strings end up by the same pattern a\a2---a k a k a k -\...a2a\ai.... 
Pictorial interpretation of the commutation relations written above one can find in [1]. Lets note 
that not all of the relations written above are independent (for details see [2], [1]). There are 
some consequences of these relations necessary, in particular, if one wishes to check validity of 
Yang-Baxter equations for the baxterization solutions, but we will not stop on these technical 
details here. 

In [1] baxterization of this algebra was successfully done. Namely, all possible solutions of the 
Yang-Baxter equation 

W i (x)W i+1 (xy)W i (y) = W i+l (y)Wi{xy)W i+1 (x) (5) 

were found in the form 

k 

Wi{x) = h + Y, ^m(x)U^ m) (6) 

m=l 

satisfying 

Wi(l) = U (7) 

in the class of rational polynomial functions a m {x). It was done by substitution of © into ©, 
which, using the Fuss-Catalan algebra, leads to a complicated system of functional equations on 
the coefficients a m (x). General solution of this system of functional equations was found to be: 

a m (x) = —x ri+r2+ - +rm - 1 (x rm -l),m = l,2,...k-l. (8) 

Pm 

ak f x ) = _^_ x n+r 2+ ...+r fc _ 2 +lZ £ (9) 

Pk-1 {I -x 

where p m = Pi(m) are defined in ([2]), r m = {±1}, and r m and other parameters entering the last 
equations, are restricted and defined by the following set of relations: r% = r^-i = 1 and either 

r q -i =r q = r k _ q = r k ^ q+l , or r g _i = -r q = -r k _ q = r k - q +i, and 

M (r s _ 1+ r,)/2 = a 2^ = a 2_ ^ ^ = ak+l q ( 1Q ) 

for any q = 2,...,k — 1. In [1] the special solution with r\ = r2 = ... = r k -\ = 1 was called 
fundamental. With this set of r q one has p = a\ — 1 = a| = ... = These equations 

were solved in the assumption that all a's are positive. We will relax this condition and permit 
two possible branches of the square root when solving (fT0|) with respect to a q . For real a q it 
means the option of some negative q's won't be ignored. We introduce the set of signs of a q , 
s q = {±l},q = 2, ...,k — 1 such that s q = s k+ i^ q . Then we get the solutions in the following 
general form 



a q (x) = \\\s 3 \ (^=) ^ J { -^^,q = 2,...,k-l, (11) 




q-l 



Er 1 ' 



= I I I4 I (z^'^zlL ( 12 ) 
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Note that an important artefact of the baxterization procedure in the case of the Fuss-Catalan 
algebras is a necessary additional requirement that all the parameters a q of different copies of 
TL subalgebras should be equal (up to some signs). Here and below we denote a\ as a - the 
only free parameter remaining after baxterization, which, we in principle permit to be negative. 
This possibility of negative a(s) will extend variety of possible physical regimes of our statistical 
models. 

In [1] the solutions with s q = r q = 1 were called fundamental. All others differ from it 
either by less powers of x appearing in a q , or alternating signs of a q . Both type of changes 
compared to the fundamental solution should respect the symmetries of parameters r and (jlOp . 
The fundamental solution can be written as 



mm = 1 - + E T^gr^r^V' + r ^—^ J l 1 « °?" (13) 



All the solutions for Boltzmann weights © we got have two important properties: unitarity 

W l {x)W l (l/x) = h (14) 



and crossing symmetry: 



Ek-2 9 
, 3=1 T 3 X 



X 



Wiixi/x) = - J ' -f—. -Wiix) (15) 
x J x(x — 1) 



X*\ fc - 2 x\ — X 



which becomes 

Wi(x 2 jx) = (-) 7* ; Wi(x) (i6) 

\ X J x[x — 1) 

for the fundamental models. Here = y/a 2 — 1, the "bar" operation means C/^ = uf" q \q = 
1, k, (uj~ = lj). In pictorial representation "bar" means rotation of TL generators by 90 
degrees, as in the crossing relation of Boltzmann weights of a lattice models. These two symme- 
tries ()14p . (|15p are very important and play crucial role in derivation of inversion relation. Here 
and below we mainly stop on few representative examples of the plenty solutions described by 
pd|) . (fr2]) . (ri3|) . We will analyze in details k = 2,3,4, considering possible regimes of the models. 



2.2 Examples and regimes 

Here we write down explicitly all the solutions found for k = 2,3,4. Two color case has only 
fundamental solution 

W?\x) = l i+ X -^U^ + 4^Lu^ (17) 
1 a a 1 — 1 — x 

In the case of three colors there are two possibilities 

Wt\x) = l t + ^±UV> ± ^2LuP ± - X2{X ~ 1) t/f (18) 

ot ay a 2 — 1 \J o? — 1 (a 2 — 1 — x) 

where the case of upper signs (3+) corresponds to what we called the fundamental solution. 
There are more possibilities in the case of four colors. It is the minimal number of colors when 
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what we called an excited solution (some r q = —1) exists. Two solutions are 



+ 



—u {1) ± x<yX - ^ 

a 1 aVa 2 - 1 
x 3 (x — 1) 



\ x ~ l ) rr(3) 



a (a 2 



(a 2 - 1) (a 2 - 1 - x) 



[7 



(19) 



The fundamental solution is (4+) and (4—) is an excited one. There are more four color excited 
solutions: 



W, 



(4'±) 



(x) 



U + 



-UP ± 



1. 



K (2) + 



o 



o 



+ 



x(x — 1) 



(a 2 



1 



-i7 



(4) 



(20) 



If we now impose positivity of all Boltzmann weights, we can find different regimes of the 
models. It is convenient to give names to the following regions on the plain (a, x): 

(21) 



I : 


1 < x < a 2 


-1,V2 < a, 


II : 


a 2 — 1 < x 


< 1, -V2 < a < -1 


III : 


< x < 1, • 


-1 < a < 0, 


IV : 


— oo < x < 


— 1, — oo < a < —1, 


V : 


— oo < x < a — 1, — 1 < a < 0, 



We will use the following parameterization of spectral and crossing parameters in the regions 
I-V: 



I, II : 


x = e u , a 2 - 


1 




III : 


x = e u , a 2 - 


1 




IV : 


x = -e u , a 2 




1 


V : 


x = -e u , a 2 




1 



e\ J : < u < A, II : A < u < 0. 



e\ A < 0, u < 
1 = e\ < u 

-e A , A < 0, A < u 



(22) 
(23) 
(24) 
(25) 



Models "2" and "4'+" are physical (have positive Boltzmann weights) in all the regimes I — V. 
"3+" is physical in the regimes I and II, and "3—" - in the regimes IV, V. "4+" has all 
Boltzmann weights positive only in the regime J, whereas models "4—" and "4'—" are never 
physical in the sense of all Boltzmann weights positivity. 

We see a rich spectrum of possibilities for different regions of parameters defining a positive 
Boltzmann weights. Only detailed investigation of the models can answer on the question what 
is the phase structure of the models and its relation to the regimes. 



2.3 Vertex and IRF representations. 

In principle it would be plausible to proceed with exact solution of these new integrable lattice 
models directly in terms of Fuss-Catalan generators. Such possibility is shortly discussed in 
the Sect. 5. Here and below we formulate integrable lattice models in terms of more concrete 
representations of the Fuss-Catalan algebra. In this context the useful fact [2] is that the Fuss- 
Catalan algebra ([I]),©,® is a subalgebra of tensor product of few copies of TL algebras with 
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different parameters TL(a.\) ® ... (g) TL(a>k), and generators E/^ may be expressed in terms of 
u (q) g TL{a q ) as 



U 



(m) 



l®...l®uf q+1) ®...®uf ] i 



® ... ® 



1 



(8) 1 i 



ode? 
even 



(26) 



Recall that one of the requirements of the baxterization procedure was the same value of \a q \ for 
different q. 

The most studied lattice models are vertex models and interaction round the face (IRF) 
models. In the first case the TL operator U- 9 ' acts on linear spaces attached to four links of 
square lattice connected to each site. Using a basis {e^} of the linear space the vertex Boltzmann 
weights W is defined as 

W(x)e<8>e = Y^Wj!f(x)ei®e j (27) 

i,j 

The fundamental representation for the usual TL generator u is known to be defined on a tensor 
product of two dimensional linear spaces in a matrix form: 



/ 1 \ 

10 

10 

\ 1 / 



it 



/ \ 

z 1 

1 1/2 

\ o o o o / 



(28) 



where z is related to the TL algebra parameter a = z + 1/z. In view of (I26p it is natural to 
define representation for the multicolored generator uj~ q ' as the tensor product of matrices (|28|) 
(with, in principle, different parameters z) and identity matrices, as in (j26|) . The linear spaces 
living on the links of the lattice are of dimension 2 fc , where k is the number of colors. 

Another representation of algebraically formulated lattice model Boltzmann weights Wi{x) 
is an IRF representation. In this case some (integer) numbers are living at each site of the 
square lattice with a constraint that on two admissible (connected by a link) sites may appear 
only numbers permitted by some fixed admissibility graph. Then the representation is given 
(see e.g. [11]) on an elementary face of the lattice with site values (i,j,k,l) in the anti clockwise 
direction 

life 



1 



'.i 



hie, u' 



(29) 



where S 1 , are components of the eigenvector of the admissibility matrix with the largest eigenvalue 

equal to the TL parameter a. a = 2 cos — , Si = sin — , with m-Coxeter number for models with 

m m ^ ^ 

ADE-Wke admissibility graphs, and a = 2 for the extended ADE admissibility graphs. Further 
construction of representation of generators goes by tensoring (I29p according to (|26j) . This 
tensor product of elementary IRF representations of TL operators will be defined on a tensor 
product admissibility graph. As we said above, not all admissibility graphs may participate in this 
tensor product, since their matrix eigenvalues are not arbitrary, but restricted by condition (jlOp . 
Possible partners among ADE models which can be tensored, were classified in [1]. Actually there 
is even more freedom - one can make a tensor product of vertex and IRF models, consistently 
fitting their parameters a. But in what follows we will concentrate on two simplest examples - 
two color vertex model and two color A n IRF models. 

The simplest possibility for vertex realization is the model "2" . In has 32 non zero Boltzmann 
weights and list of them one can find in [1]. Looking at this list, one can see that the model 
is not of ice type - there is no a natural "charge" conservation in each vertex. Linear space 
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living on each of four edges of a vertex is four dimensional, such that R-matrix is 16 x 16. An 
interesting observation is that it is possible to attach "charge" to these four spaces in a way 
preserving "charge" conservation, but on a cost of charge degeneracy of the spaces. For instance, 
the attachment: space — ► charge 0, space 1 — ► charge 1, space 2 — > charge 1, space 3 — > charge 
3 - leads to "charge" conservation in each non zero vertex. 

We formulate the model in the language of spin chain models, and consider the Boltzmann 
weights W as an L-operator. (Recall that R-matrix is related to L-operator by permutation P 
of outgoing spaces: R = PL, or Rgjj = Wxh- We define the transfer matrix as 



T, 



0( e)^) = (w 0{e) (x))zz (WW*)):: (w oie) ( X )) a c ^ ... (w eio) ( X )) 



OJV-lCJV 
CN-lbN 



(30) 



There is a summation over repeated indices here, and indices o/e mean oddness (evenness) of the 
south site along the line of the auxiliary space c... As we see, T is an alternating product of 
Boltzmann weights on even and odd sites. (Recall that in the multicolored TL models they are 
different). In the language of spin chains Boltzmann weights may be considered as L-operators. 
Unitarity condition looks like 



5 a e 8) 



(31) 



There is an important property of BW, which is a mixture of crossing (I15p . (|16p and unitarity 
relation: 



W, 



o(e)l 



(32) 



Here t% means transpose in the second space of Boltzmann weights, corresponding to the auxiliary 
space of the transfer matrix. The matrix acts in the first space and has the form 



M 



( o 



V iA 2 






1/z 




1\ 

1/z 



J 



(33) 



Another important symmetry is 

[W o(e) (x)) 



ab 
cd 



PW e(o) (x)P = (W e(o) (x)) 



6a 
dc 



(34) 



The relations (|3ip . (l32p are almost enough for derivation of the inversion relation (see below). 

As the simplest example of two color IRF models one can consider the case of tensor product 
admissibility graph A n x A n . The general form of the RSOS Boltzmann weights are defined by 

021,(21}, 



W e (o)(ii' ,jf ,kk' \ll'\x) = 5 ik 5i'k' + 



a 



fii'k'8jl^-g 



+ 



x(x — 1) 



1 



OikOj'V — s~, — 

SjSf 



even site 
odd site 



+ 



(35) 



Pairs of variables (£,£') , (j, f) , (k,k') , (1,1'), siting at the south, east, north and west corners 
of a face, are any admissible set of integers in the range from 1 to n. Each of two quartets 
(i,j,k,l) and (i' , f , k' , I') are any permitted by A n admissibility graph diagram independently, 
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• — * 



Figure 1: (A n ) 2 admissibility graph decomposition: a (upper) n is even, b (lower) n is odd 

and Sk = sin , a = 2 cos • One can explicitly check that these Boltzmann weights satisfy 
the following unitarity and crossing relations 

^2 W e(o) (ii',jj',kk',ll'\x) W o(e ) (kk',jj',mm',ll'\l/x) = 5 im 5 ilm , (36) 

k,k' 



W e(o) (ii',jf,kk',ll'\xl/x) = a \, X \ S Jc lS J ,S J' W o{e) (jj',kk',ll',ii'\x) (37) 

X[X — I) V DiDkOi' >Jk' 

An interesting and important feature of the Boltzmann weights (|35p is their admissibility 
graph decomposition property. One can check that admissibility rules coming from (|35p lead to 
two independently existing subsets of site variables, each with its own admissibility graph. The 
form of these two graphs depends on oddness of n (see Fig. la, Fig. lb). 

Again, the basic object of investigation will be the transfer matrix 

N/2 

T o(e)(x)ffJ } = JJ H / " o(e) (a i a-,aj + ia- +1 ,/5i + i/3- +1 ,/5i/5-|j;) x 

i=l 

W e ^(a i+1 a' i+1 , a i+2 a' i+2 , (3 i+2 (3' i+2 , (3 i+ i(3' i+1 \x) (38) 

where the trace identification ao = «o = a Ar>A) = PniP'o = f3' N is imposed. Recall here an 
important and well known difference in definition of TM for vertex and IRF cases. In the vertex 
case the auxiliary space of TM contains a sum over thermodynamically many (A^) variables Cj in 
(|30p . whereas the IRF TM does not contain any auxiliary space. 
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3 



I 



4 



6 3 2 
• •- 



10 



9 



16 



a. n=3 



Figure 2: Admissibility graph decomposition for (.A3) 2 and {A^Y 



Let us look in more details at the two color (A3) 2 graph IRF model. This model, in a sense, 
may be called A3" two color" Ising model, since, as it is well known (see e.g. [12]), A3 graph IRF 
model is equivalent to the Ising model. Recall that A3 admissibility rule requires that the height 
value is fixed and equal to 2 (the value of the middle site of A3 graph) on the even sublattice 
of two dimensional square lattice, whereas the heights on the odd sublattice take the values 1 
and 3, forming an Ising model. The situation is different in the two color case and the model 
seems to be more complicated then Ising model. As one can see by explicit analysis of non zero 
Boltzmann weights for the case n = 3, the RSOS model splits into two subsets of variables living 
on the sites of a face. If one denotes ({a}) the pairs of two numbers at the sites of a face as 

(ij)=a = 3(i-l)+j (39) 

the decomposition of the admissibility graph into two ones looks like it is shown on the Fig. 2. It 
means that there are either even or odd heights in the new notation, living on each site of the 
model. It means the model factorizes into two separate submodels. 

In the same way, admissibility graph of the two color (A4) 2 graph IRF model, after the 
recounting of the site variables 

(ij) = a = 4(i-l)+j (40) 

decomposes into two disconnected graphs (see Fig. 4a, b), dividing the model into two submodels - 
one with the site variables a = {1, 3, 6, 8, 9, 11, 14, 16}, another - with a = {2, 4, 5, 7, 10, 12, 13, 15}. 
Again there are two separate submodels. 

Note that the relations ([3?]) are valid for each of two submodels of (A n ) 2 admissibility graph 
decomposition. 

3 Inversion functional relations and their solutions. 

The most full and contemporary investigation of the models described in the previous section 
should probably follow the standard methods of algebraic Bethe ansatz. Unfortunately the ob- 
structions we met on this way didn't allow us to proceed in this direction, and algebraic Bethe 
ansatz for these models remains an open problem. But some information about ground state 
energy and free energy, spectrum of low lying excitations, and even about some correlation 
functions, can be extracted using a method known for a long time [7], which doesn't require 
a knowledge of explicit form BA equations and transfer matrix eigenvalues. It is a method of 
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inversion functional relations on the transfer matrix eigenvalues, which was successfully applied 
to few integrable lattice models [10]. It would be interesting to develop a universal algebraic 
method of derivation of inversion relations in terms of multicolored TL algebra generators, as it 
was partly done in [13] for usual TL algebra models, but in what follows we are going to consider 
the inversion relation in vertex and IRF representations of multicolored TL algebra, which we 
discussed in the previous section. 

3.1 Vertex model 

We start from the vertex model. Consider the product 

(Mf N T (x*x)(Mf N T e (x) (41) 

Since transfer matrices commute, their eigenvalues will satisfy the same functional relations as 
transfer matrices themselves. The above expression explicitly looks as 



Tr A Tr B 



mCi W Ia Ci (xlx)M Cl M c * W\ A AC2 (xlx)M c * ..a 



W eBCl (x)W oB c 2 (x)...] 

Insertion of the equal to 1 factor {5 ai b t + 1 — ^fej i n any place of this product leads to a split 
of the above expression into two terms - the first corresponding to the <5 a A> an d the second - 
to (1 — Sa-fy). Due to the relations (|32p . (|3ip the first term turns out to be proportional to the 
identity operator in all the Cj spaces. Oppositely, the second term is not giving an identity 
operator in any of C% spaces. According to (|32p . the proportionality coefficient for the first term 



will be -^w ( x-^x-i-i) J • O ne can snow that the expansion of the second term in powers of x 

starts from the x N ^ 2 , therefore in the region where x > 1, e.g. in the regime I, the second term 
is exponentially small compared to the first term in the thermodynamic limit. This leads to the 
following inversion relation 

(Mf N T (xlx) (M f N T e (x) = w "-?_!) ) J Ci ® J Q» ® Ic " + °( e ~") < 42 ) 

One can also see that specific of the matrix M is such that the relation between highest eigenvalue 
of {M) m T {xlx) (M)® N and of T (x 2 x) is by factor It leads to the following functional 

relation for eigenvalues of transfer matrix in the thermodynamic limit N — > do, when the terms 
0{e~ N ) can be neglected 

, T 2 „-i \N 

2 \ k / \ I X^ 



A (^) A w = U'V--i) J (43) 

Another relation is the unitarity relation. If one considers the product 

T (x)T e (l/x) = Tr A Tr B [W oA c 1 ( x )WeAC 2 (x)...W eB c 1 (l/x)W oB c 2 (l/x)...} (44) 

then by insertion of the space permutation operator P (P 2 = 1) between the Boltzmann weights 
of, say, the second TM, and using (f34"|) . one gets 



E ^o(x)**W e (x)J«...W' (l/x)^W'e(l/x)5;g... (45) 

{a},{b},{c} 



10 



If one inserts one delta symbol S ei+iai for some i, then due to the relation (f3Tj) all the product 
will be reduces to the identity operator Id = 6 ei d 1 ---S eN d N - One can also show that other terms 
successfully cancel. This leads to the second functional relation on the transfer matrix eigenvalues 

A(x)A(l/x) = 1 (46) 

The functional equations (|43|) . (|46j) can be solved by standard methods. These equations in 
the parameterization (122j) for regime /, can be written in terms of A(u) = A 1 ^ N (x), as 

A(u)A(\ + u) = (47) 
A(u)A(-u) = 1 

In order to solve uniquely this system of functional equations, one needs additional information 
about analyticity of the function A(u). One can look for solutions in the form A(u) = F(u)f(u), 
where F(u) is a general analytic non zero (ANZ) in the quadrant < Re(u) < A, < Im(u) < 2ir 
solution of the system 

F(u)F(X + u) = (48) 



and f(u) is a solution of 



F(u)F(-u) 



/(«)/(« + A) = 1 (49) 
f(u)f(-u) = 1 



with a given set of poles u p and zeros u z in the same quadrant. We have cut the strip < 
Im(u) < 2ir in the imaginary direction because of the obvious 2iri periodicity of solutions of the 
system (j47]). The cut in the real axis direction < Re(u) < A is a consequence of the relation 
f{u + 2A) = f(u), which one can easily derive from (|49p . 
The general ANZ solution of (J48]) has the form 



oo ( (2j+l)A+« _ -A C (2j+2)A-u _ l) 

f <"> = - » n \ e2J > + ,_ 1)( lU»A-.-i ) ) < 5 °) 

For the double periodic function f(u) we have the system of functional equations (|49p . Pro- 
vided we know the set of its poles and zeros in the periodicity quadrant, the solution is uniquely 
fixed by Liouville's theorem: 

f(u) = ± PJ ^/K^snh (^^-(u - u z )^j Y\ y/Kisnh {^^-{u -u p - X)j (51) 

where snh is a standard elliptic function with modulus K\ defined by the requirement that 
corresponding quarter periods K\ , K[ are related by = ^ . (j50|) and (f5TI) give general form 
of transfer matrix highest eigenvalue as a function of spectral parameter u, provided one knows 
positions of its poles and zeros in the periodicity region. Obviously transfer matrix eigenvalues 
do not have poles in the physical region. Some numerics we did for the lattices of small size in 
regime / support the conjecture that the ground state eigenvalue do not have also zeros in the 
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periodicity region, and f(u) = 1 for the ground state. Using (|50p one can calculate the ground 
state energy of the model 



F ~ d 
au 



(lnA(u)) n = -Y : (52) 

V V JJu=0 ^sinh(j + l)Asinh(j + i)A V ; 



3.2 IRF model 



Discussing the family of IRF models, we will concentrate on two color (A n ) 2 models. In the same 
way, the product of two transfer matrices (|38[) for the IRF model has the form 

N/2 

T (xlx)T e (x) = ^2 n w o(ai,ai+i,Pi+i,Pi\xlx)We{a i+ i,a i+ 2,f3 i+ 2,Pi + i\xlx) * 

{oti} 8=1 

l^e(7i>7i+l>^+l ) «i|«)^o(7i+i ! 7i+2>«i+2,ai+l|ic) 



Using crossing and unitarity properties f|3T[) of the BW one can easily see that insertion of one 
Kronecker Sarf- f° r some one i into the last expression, reduces it to the identity operator for all 



„2_ ,.-1 



N 



i with the prefactor y x -i^ x -i_^ J • R means that the identity element is in a sense "decoupled" 
from the others. With some effort, one can show that, like in the vertex case, x expansion of 
operators different from identity starts from x in the maximal power less then N. It leads to the 
same functional relation as (|43p . 

2 -1 \ N 
X* — X \ „ , 



T a (xix)T e (x) = ^ x _* _ \ +0(e^) (53) 

Unfortunately in the IRF case the situation is different with the unitarity relation (|46p . Were we 
try to generate an identity operator for product of two transfer matrices T (x)T e (l / x) , we would 
need to introduce 2N delta symbols in order to make the relation (|36|) working, without involving 
the crossing relation (|37p . (Recall that there is no any auxiliary space trace in the definition of 
the IRF transfer matrix). Apparently such drastic deformation operator cannot be negligibly 
small in the thermodynamic limit N — > oo. There is the well known way to avoid this obstacle: 
one can pass to the bigger object - partition function. 

Consider partition function for 2N x 2M square lattice. It can be written in different ways 

Z(x) = £ (T^(x)T^(x)) 2M = £ (T^(x)T^{x)) 2N 

p.b.c. p.b.c. 

where the transfer matrix product means summation over the repeating site variables of the 
multiplied transfer matrices, and ^ P f> c summation over periodic boundary conditions in both 
directions. One can introduce a set gM+N of delta symbols matching the boundary configurations 
for both statistical sums (see the dashed lines on the illustrating Fig. 3 left). Insertion of the 
identity operator in the form 1 = gM+N + (1 — gM+N) separates the partition functions product 
into two pieces: 

2N / ,«« \2M 



Z(l/x)Z(x) = £ (T^\l/x)Ti M \l/x)) gM+N {T^ N \x)T^\x)) + 0(x) (54) 

One can see that the first term "collapses" into operator identifying configurations on all the sites 
of both lattices due to the relation (|36|). This "collapse" goes plaquette by plaquette, without 
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M 



2 



M m y=x X 

Figure 3: Illustration of inversion relation derivation for IRF model 

any relation to the transfer matrices. The second term 0(x) represents deviations from identity. 
Here we will assume that the boundary changing operator qm+n produces only subdominant 
terms for partition function. This assumption was shown to be correct for some IRF models [10], 
but investigation of this question is needed in the case of our IRF models. With this assumption 
the first term in (|54p is dominating over 0{x). 

In the similar way one can translate the relation (|53h valid by itself, to the language of the 
partition functions product. One can see (for illustration see the Fig. 3 right) that the same gM+N 
does the job: 

Z(xlx)Z{x) = {T^\xlx)T^ N \xlx)) 2M g M+N (^(^(^ + 0{x) (55) 

In the first term, as it is illustrated on the Fig. 3b, the product of the right column transfer 
matrix of Z{x\x) with the lowest row transfer matrix of Z[x) by (|53p is proportional to the 
identity operator. This "collapse" to the identity continues to the product of the next to right 
column with the next to lowest row, and so on. By the same reasons as in the previous (|54p 
relation, one can believe that the first term is dominating in the thermodynamic limit. 

The dominating terms of equations (|54p . (|55p lead to the system of functional equations for 
partition function density z(x) = Z(x) 1 ' MN 

z(x)z(l/x) = 1 (56) 

which coincides with the equations (|47|) in the proper regime parameterization. Numerical sim- 
ulations show that for two color (A3) 2 and (A4) 2 models z(u) is ANZ in the periodicity domain. 
It means that the solution for partition function in the ground state z(u) is given by (|50p . An 
interesting property of it is that z(0) = 1. 
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3.3 Spectrum of low lying excitations 

In the same way one can consider the low lying excitations, i.e. those states for which the transfer 
matrix eigenvalue A(u) is finitely different from the ground state one A gs (u) in the thermodynamic 
limit. For instance, define finite function 

l(u) = lim -^L (57) 

for vertex model. Again l(u) can be represented as a product of ANZ function F(u) and another 
one f(u) satisfying the same functional equations (j4*8|) . (|4l)|) . But this time f(u) / 1, and some 
zero(s) u z enter into the periodicity domain. The excitation spectrum can be calculated from the 
relation between energy and momentum of excitations 

E-Eo = -^^(lnV^)) u =o = £e(p(u*)) 

{Uz} {Uz} 

P _P = -f^lnV(O) = ^p(n,) 

{u z } {Uz} 

where lmp(u) is the contribution of each snh factor in the product (|5ip . Using (|5ip we get the 
dispersion relation: 

e(p) = K'JJl - k) 2 + AKsm 2 p (58) 

with non zero mass gap 

A = K'(l - k) (59) 

So, at least in the regime /, where one can neglect the additional terms in the inversion relation, 
all the excitations are with a non zero energy gap. It means there is no conformal critical point 
in this regime. 

Importance of the regime in all what was done above, is crucial in two aspects. First, in order 
to neglect the non identity terms in the inversion relation one needs |x| > 1, which is correct not 
only for regime /, but also for regime IV and x < — 1 subregion of the regime V (we will call 
it V). Second, the numerics, which is necessary to be sure in ANZ condition, essentially differ 
for different regimes. The analog of (I47j) for vertex model in a parameterization (I24l ). (|25j) of the 
regimes IV, V looks like 

A(u)A(A + u) = e -^f (60) 
A(u)A(-u) = 1 

where in the regime IV p = — 1, and in the regime V p = 1. Equation (j49[) remains as it is. It 
means the solution for F in this case is given by 

oo ( (2j+l)X+u _ \ ( (2j+2)X-u , i\ 

^-'■■ti) ^,.,,^^---,) (61) 

f(u) is given by the same solution (|5ip . With the assumption of ANZ condition for A(n) it leads 
to the following ground state energies of the two color vertex model in the regimes IV and V: 



4 v = -y^ 



cosh l E v = _ cosh l m 

cosh(j + l)Acosh(j + |)A' cosh(j + l)Asinh(j + \)\ 



14 



The analysis of low lying excitations remains the same, giving the same finite energy gap (|59D . 

In the same way one can modify the solutions for partition function z(u) of two color IRF 
models in regimes IV and V' . With an ANZ hypothesis, it is given by the same expression (|6ip 
with the same property z(0) = 1. 



4 Other models. 



A natural conjecture is that all what was done in the previous section for two colored models, 
can be generalized for what was called above fundamental models with any number of colors in 
the regime analogous to regime /. Lets point out that both three and four colored models have 
only regimes I and II as regimes with all Boltzmann weights positive. As we said, x > 1 only in 
the regime /, which was crucial for derivation of inversion relation. Instead of (|32p we have 



Li(x) 



fc-2 X (X - 1) 



X 



(63) 



with some matrix M acting in the first space. Together with (|3ip it gives the following inversion 
relations 



T(x)T(xlx) = {x*x) k ~ 2 
T(x)T{l/x) = Id 



x-^x- 1 - 1) 



N 



ld + 0(e 



-N\ 



As we see, they differ from (|53|) is minimal way - by a prefactor. The same procedure as previously 
gives 



A(«) 
A(«) 



A(«)/(«) 



e (k-2)u ,u 



3 (2j+l)A+« 



> 11 ( e 2jX+u _ 1) ( e (2j+l)A-« _ ^ 



with the same expression (|5ip for /(«). It leads to integer shift of the ground state energy in 
k — 2, compared to (f52l) . provided still / = 1 (ANZ hypothesis). This of course should be checked 
numerically. 

Looking at the other models which have some physical regimes, one can believe that the 
general form of inversion relations both for fundamental models with some negative signs and for 
excited models with or without negative signs, is just 



T{x)T{xlx) 
T[x)T(l/x) 



Id 



X 



n N 



x-^x- 1 - 1) 



ld + 0(e 



Form of solution of these equations depends on regime. We see that different models have the 
same set of functional inversion relations, and sometimes even the same set of regimes. 

One can believe that if there is a difference in ground state or low lying excited states between 
these models, then it is expressed in the amount and location of poles and zeros of function f(u). 
An information about these poles and zeros may come from numerics. 
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5 Discussion 



In this paper we made an attempt to analyze vertex and IRF multicolored lattice models. We saw 
an interesting structure of admissibility graph for IRF models. Being unable to solve the models 
by algebraic Bethe ansatz, we tried to analyze them using the inversion relations. The main result 
of this attempt - all the models have a mass gap in the considered regime. Another evidence for 
absence of conformal points in this regime we saw formally trying to construct new link invariants, 
following a standard procedure of extracting a braid group representation from the Boltzmann 
weights either in vertex or in IRF form (see [14] and references therein). The objection we met 
in this way is that it is impossible to get a well defined braid group representation taking a limit 
u — ► ioo, or some other one, for the two color Boltzmann weights W(u). Also the formal sufficient 
condition for construction of link invariants from W(u) one can find in [14] is not satisfied. But 
one cannot reject a possibility of gapless behavior for the models we considered in other regimes, 
as it was pointed out in [6] for dimerised coupled spin chain. 

Let us stress an almost unrestricted amount of possibilities in construction of new integrable 
multicolored TL models demonstrated above. Certainly these zoo of new models deserves a 
classification which, as far as we know, is missing for today. 

An interesting possible alternative is coordinate Bethe ansatz formulated in terms of the Fuss- 
Catalan algebra itself. Such kind of Bethe ansatz formulation, purely in terms of usual TL algebra 
generators for the XXZ model, was applied in [15]. The basis for wave function representation 
was chosen as a basis of specially defined two sided ideals of words of TL generators. It would 
be interesting to apply such formulation to the multicolored TL lattice models. 

This paper is only a first touch to these new integrable models. Development and applica- 
tions of known methods of integrable models, such as algebraical Bethe ansatz, transfer matrix 
functional relations leading to a (system of) integral equations, a role of Q-operator in these 
models, and so on - are open and interesting problems. 

6 Acknowledgements 

A.B. is grateful to M.J.Martins for very useful communications. We are thankful to Einstein 
Center for financial support. 

References 

[1] P. Di Francesco, Nucl. Phys., B532, (1998), p.609, [hep-th/9807074] ; 
[2] D. Bisch and V.Jones, Inv. Math., 128, (1997), p.89. 

[3] P. Di Francesco, E.Guitter, Phys. Rep., 415, (2005), p.l. jcond-mat/0505293] . 

[4] M.J.Martins, "Lorentz Lattice Gases and Spin Chains", Proceedings of the APCTP Winter 
School, Cheju Island, Korea, Eds. Changrim Ahn, Chaicho Rim and Ryu Sasaki, (2000), 
p.179. 

[5] M.J.Martins, B. Nienhuis, J. of Phys. A -Math. Gen. 31 (1998), p. L723. 
[6] M.J.Martins, B. Nienhuis, Phys. Rev. Let, 23 (2000), p. 4956. 



16 



[7] R.Baxter, J. Stat. Phys., 28, (1982), p.l. 

[8] R.Baxter, "Exactly solved models in statistical mechanics" , Acad. Press, London, 1982. 

[9] RPearce, Phys. Rev. Lett, 58, (1987), p. 1502. 

[10] A.Kluemper, J.Zittartz, Z. fur Phys., B71, (1988), p.495; 
A.Kluemper, Europhys. Lett, 9 (8), (1989), p.815.; 
A.Kluemper, J. of Phys., A23, (1990), p.809. 

A.Kluemper, A.Schadschneider and J.Zittartz, Z. fur Phys., B76, (1989), p. 247; 
[11] V. Pasquier, Nucl. Phys., B285, (1987), p.162. 
[12] P. Di Francesco, Int. J of Mod. Phys., A7, (1992), p.407. 

[13] P.P.Martin, "Potts models and related problems in statistical mechanics" , Singapore: World 
Scientific, 1991. 

[14] T.Deguchi, M.Wadati, Y.Akutsu, Adv. Stud, in Pure Math., 19, (1989), p.193. 

[15] D.Levy, Phys. Rev. Lett, 64, (1990) p.499; 

P.P.Martin, H.Saleur, Int. J. of Mod. Phys., B8, (1994) p.3637. 



17 



